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ABSTRACT 

We present narrow and broad K-band observations of the Class O/I source 
IRAS 18148-0440 that span 17 years. The infrared nebula associated with this 
protostar in the L483 dark cloud is both morphologically and photometrically 
variable on a time scale of only a few months. This nebula appears to be an 
infrared analogue to other well-known optically visible variable nebulae associated 
with young stars, such as Bubble's Variable Nebula. Along with Cepheus A, 
this is one of the first large variable nebulae to be found that is only visible 
in the infrared. The variability of this nebula is most likely due to changing 
illumination of the cloud rather than any motion of the structure in the nebula. 
Both morphological and photometric changes are observed on a time scale only a 
few times longer than the light crossing time of the nebula, suggesting very rapid 
intrinsic changes in the illumination of the nebula. Our narrow-band observations 
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also found that H2 knots are found nearly twice as far to the east of the source 
as to its west, and that H2 emission extends farther east of the source than the 
previously known CO outflow. 

Subject headings: ISM:individual(L483), ISM:reflection nebulae, ISM: structure, 
ISM: kinematics and dynamics 



1. Introduction 



Variable nebulae, well-known examples being Hind's Variable Nebula, R CrA, Bubble's 
Variable Nebula, and PV Cephei, are now known to be associated with embedded young stars 
and were discovered through visible light observations. Hind's Variable Nebula (Hind 1864), 
disco vered near T Ta uri in 1852, was easily visible in 1855-56, but was not seen by Struve in 
1861 jSarnardl Il895l ). However, Struve observed a small nebula 4' to the west that had not 
been previously v isible. The nebula around R CrA was found to be variable before 1890. 
Knox Shawl Jl916h describes the changes in the morphology of the nebula, how the changes 
in the brightness of the nebula are uncorrelated with the brightness of R CrA itself, and that 
change s are not i ceabl e on a time scale of a week but are not noticeable between consecutive 
nights. iHubbld (119161 ) found that several features in NGC 2261 ( Hubble's V a riable Nebula) 
were observed to vanish, then later reappear in the same location. iLamplandl (119261 ) success- 
fully explained the changes in morphology as changes in the illumination of the cloud rather 
than on account of any actual motion of the structure in the nebula. PV Cephei is another 
young star with an optically visible variable nebula (Cohen et al. 1977; Gyul'budagyan & 
Magakyan, 1977). Morphological changes are seen on photographic plates within a month, 
and the point-like sou rce at the southern tip of th e nebula varied by over four magnitudes in 
less than six months. iHodapp &: Bressert I (120081 ) found that the infrared reflection nebula 
illuminated by Cep A also shows photometric and morphological variability. Changes in the 
illumination of the circumstellar m aterial were proposed to account for the observed changes 
in morphology. McNeil's Nebula (iReipurth &: AspinI |2004J ) is a recent example of a large, 
optically visible variable reflection nebula that became visible when V1647 Ori experienced 
an eruption. Optical photometry (IBriceiio et al.l |2004J ) shows that the photometry of the 
nebula closely followed th e photometry of the central star, and that the morphology of the 
nebula (IFedele et al.l 120071 ) was largely static during the eruption. 



IRAS 18148-0440 was identified as a candidate protostar by iParkerl (119881 ) on account of 
its IRAS colors and coincidence with the L483 dark cloud. This source has been classified as 



a Class YSO bv lAndre et all mm . and as a Class I YS O bvlPezzuto et all (120021 ) based 
on its ISO 60 jim — 100 /xm vs. 100 /im — 170 /im colors. iTafalla et al.l (I2OOOI ) suggest that 
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this source is in transition between Class and I on account that it has an SED characteristic 
of a Class source but its outflow properties and near-IR reflection nebula suggest that it is 
a Class I source. Compared to Cep A, the source in L483 is a lower mass, lower luminosity 
object that is likely to be at an earlier evolut ionar y state. The distance t o L483 has been 
given as 2 pc by iDame fc Thaddeud (119851 ) and iGoodman et al.l (Il993l ). and 250 pc by 
Felli et all Jl992h . 



Our observations have found that IRAS 18148-0440 is the source of a new variable nebula 



that is only visible in the infrared. [Fuller et al.l (119951 ) presented images of the infrared nebula 
at J, H, K, and H2, demonstrating the alignment of the H2 knots with the CO outflow and 
the coincidence of the near-IR reflection nebula with the blue shifted lobe of the outflow. 
They estimate that the dynami cal age of the outflo w is 1.3x10^ years, and that the outflow 
is tilted 40° to the line of sight. iBuckle et al.l (119991 ) studied the excitation temperature and 
velocity of the shocks associated with the outflow using long slit spectroscopic observations 
of the H2 emission frorn knots in the jet. Derived properties of L483 are summarized by 
Fuller fc WoottenI (120001 ). who observed submillimeter dust continuum emission and NH3 
line emission. They interpret the NH3 line widths and velocities as a signature of infall, and 
present a model to account for their observations. 

In this paper we present K-band images of IRAS 18148-0440 collected over the span of 
17 years, as well as an image that shows the full extent and asymmetry of the H2 knots. 
The observations show that the infrared nebula is both photometrically and morphologically 
variable. Our goal is to present a time series of K-band images and photometry of this nebula 
to demonstrate its photometric and morphological variability. 



2. Observations and Photometry 



We observed IRAS 18148-0440 (^18^ IT'" 29 8 " -04 ° 39' 38" J200 0) in L483 on 9 occa - 
sions, and include observations from lFuUer et al.l (Il995l ) and 2MASS (ISkrutskie et al.ll2006l ) 
for a total of 11 e pochs . The location of the centimeter source in the nebula given by 
Fuller fc WoottenI (I2OOOI ) is less than an arcsecond from the position noted above. The 



observations used in this paper are summarized in Table 1. The data acquisit i on an d re- 
duction for t he dat a taken in Oct. 1990 and June 1 993 are desc ri bed i n lHodappI (Il994j ) and 
Fuller et al. (1995), respectively. Observations by Fuller et al. ( 1995[) were t aken through 
Barr near-infrared filters. The data taken with QUIRC (iHodapp et al.lll996l ) at K and K' 
( Wainscoat fc Cowielll992l ) and the SpeX guider (IRayner et al.ll2003l ) are in the MKO pho- 
tometric system (Simons & Tokunaga 2002, Tokunaga & Simons 2002). 
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All data from QUIRC and SpeX were taken while dithering, usually with an equal 
number of dithered sky frames offset from the target. For the data reduction, a dark frame, 
of the same exposure time as the data, was subtracted from each frame. A sky flat was 
constructed from the dithered images of the sky near the target by first scaling these images 
to the same median value, then median combining them together with a min-max rejection. 
We divided each image of the target by a normalized sky fiat, then subtracted the median 
sky level from each image. The resulting images were then aligned and averaged together 
with a sigma clipping rejection. 

The observations of L483 through the H2 v=l-0 S(l) 2.121 /im filter have a slightly 
different plate scale than the K-band images. Before the K-band image could be subtracted 
from the H2 image, we used the IRAF tools geotrans and geomap to match the star fields. 
We then used a Gaussian convolution to blur the K-band image to match the FWHM of 
the H2 image, scaled the resulting image so that the stars had the same peak counts, then 
subtracted the scaled K-band image from the H2 image. 

We measured the photometry of the L483 infrared nebula relative to four nearby field 
stars since some of the nights were non-photometric. The photometry of the field stars was 
measured using an aperture 1.5 times the FWHM to ensure that the same fraction of the 
stars' light was in the aperture despite the widely varying angular resolution of our data set. 
Since the nebula is spatially extended, and elongated east to west, we used a rectangular 
aperture measuring 40" east to west and 20" north to south to measure the photometry of 
the nebula. The average of the counts from the nearby sky was subtracted from the image 
before the counts in the rectangular aperture were summed. 



Discussion 



3.1. Classification 



Class YSOs are characterized by having an SED that is a cool single temperature 
blackbody (lLada]ll99ll ). with little or no flux in the near-IR, whereas the SED of a Class I 
YSO can not be fit by a single temperature blackbody and has significant flux in the near-IR. 
As mentioned in the introduction, various authors have classified this object as a Class 
YSO, a Class I YSO, or an object in transition between the two. We combined data from 
2MASS (J, H, and K), Spitzer (3.6 /xm, 4.5 ^m, 5.8 /im, 8.0 /im), IRAS (12 /im, 25 /xm, 
60 /im, 100 fxm) and JCMT/SCUBA (850 /im. Fuller & Wooten 2000) to create the SED 
shown in Figure 1. We note that the 12 /im flux is only an upper limit. 



We chose to fit the IRAS fluxes with a blackbody function since that wavelength range 
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is least affected by ffux from the cold envelope and is least attenuated by extinction. When 
these data are fit with a 45 K blackbody (the best-fit temperature), the 850 fim flux measure- 
ment is much greater than predicted by the fit. If the 850 /xm flux measurement is iii c luded , 
the best-fit temperature is 35 K. The 850 fim map presented by iFuUer &: WoottenI (120001 ) 
shows emission extending ~60" from the source. As such, much of the emission at this wave- 
length is likely to be from the cool envelope around the source rather than the source itself. 
Despite the large amount of extinction towards this object, there is still significant amount 
of flux at wavelengths shorter than 10 ^m. We did not fit a separate blackbody function to 
the flux measurements at wavelengths shorter than 10 fim since they are strongly attenuated 
by extinction and also strongly affected by scattered light. This excess flux at wavelengths 
shorter than 10 /im supports the conclusion that L483 is not a Class YSO, but has more 
in common with deeply embedded Class I YSOs. 



3.2. Change in Brightness 

Accurate relative photometry of this object is very difficult. First, the object is large, 
diffuse, and irregularly shaped. While there is a dearth of field stars immediately around 
the nebula due to the dark cloud that embeds the source, the nearby field stars can affect 
the photometry if they lie either in the aperture or in the region used to fit the sky level. 
Second, the observations were not part of a focused effort, and were made over 17 years 
using 5 different cameras and filters from at least 3 different photometric systems. However, 
the different filters and color responses of the cameras affect our relative photometry of the 
nebula only so far as the color of the nebula is different than the color of the field stars used 
as the photometric references. The mean 2MASS H— K color of the reference stars we used 
is 1.66, which closely matches the 2MASS color of the nebula (H— K=1.85). Our H and 
K observations in July 1999 show that the nebula has a H— K color of 2.06, similar to the 
2MASS color. 

The symbols in Figure 2 represent relative photometry of the whole nebula using differ- 
ent field stars. In some cases, only one star was within the field of view that could be used 
as a photometric reference. Since the relative photometry of the nebula from different field 
stars are very close to each other, the color effects described in the previous paragraph do 
not appear to significantly affect our photometry compared to the photometric variability 
of the nebula. Furthermore, using the observed colors to transform the photometry to a 
common photometric system can be misleading since the color of the nebula may also be 



variable. iFuUer et al.l (119951 ) state that the K-band peak of the nebula has an H— K color of 



4.1, whereas the H— K color of the brightest part of the nebula in our July 1999 data is 2.30 
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in an 8" aperture. Variability of the color of the nebula may be due to a combination of 
the color gradient across the nebula (the nebula is bluer towards the west) with the changes 
in the brightness of various knots in the nebula. With our current data set, we cannot 
determine if individual knots have a variable color. 

Over the 17-year span of observations, the brightness of the nebula changed by over 
two magnitudes, and changes were observed on a time scale comparable to the light crossing 
time of the nebula. The integrated K-band magnitude of the nebula varied from 9.8 to 12.0 
(Figure 2), and changes in brightness by about 1 magnitude were observed within a two- 
month time span. It is likely that significant changes in the brightness of the nebula occur 
on much shorter time scales than the spans between our observations. We did not observe 
frequently enough to determine the minimum time scale of the variability, or to determine 
if the variability is periodic. The center of the nebula at K-band is ~ 15" away from the 
source, or 3000 AU at a distance of 200 pc. As such, the light crossing time is ~ 17 days. 
The illumination of the nebula must be significantly and rapidly changing to show significant 
brightness changes on a time scale of, at most, four times longer than the light crossing time. 
Our photometry also shows that the nebula appears to be fading at an average rate of ~ 0.08 
magnitudes per year. 



3.3. Morphological Variability 

We noticed changes in the apparent morph ology of the nebula over the 17- year span of 



our observations. Although we have images from lFuUer et al.l ( 1l995l ) from 1993, 2MASS from 
1999, and our data from 2002, these are not included in Figures 3 and 4 since they do not 
have sufficient angular resolution to show changes in morphological details. The 8 epochs of 
data that are used in the following analysis were added together to create a master image of 
the nebula that is representative of its average state. Figure 3 labels several regions that were 
observed to change. Figure 4 shows a K or K' image of the nebula as well as the difference 
between the image taken on a particular date and the master image. For each difference, the 
input image was scaled to match the brightness of the master image to compensate for the 
overall variability of the nebula and to better show the morphological changes. The angular 
resolution within these figures varies by over a factor of 2, and is documented in Table 1. 

The differenced images in the right column of Figure 4 clearly show the morphological 
changes in the nebula. In 1990, the central part of the western lobe of the nebula was 
brighter, whereas the opposite was true 9 years later. In July 1999, the southern half of the 
western lobe was unusually bright, but two months later this had changed. Region A was 
brighter than average in 2003 only. By 2004, region B was unusually bright, which was also 
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observed in 2007. The excess illumination of the nebula in 2007 appears to point away from 
the expected location of the source (roughly between regions A and C). 

The morphological features in the nebula tend to vary in brightness together, but can 
change brightness independently of each other. Light curves for the four regions outlined 
in Figure 3 are shown in Figure 5. Region A appears to brighten significantly from 1999 
to 2003. However, its integrated K-band magnitude brightens by only ~ 0.2 magnitudes 
whereas the other three regions fade by ~ 1 magnitude. Region A is typically fainter than 
the other three regions, and became much fainter than the other regions between 2004 and 
2007. Regions B and D had a similar brightness until May 2004, when regions C and D faded 
by ~ 0.6 magnitudes but region B stayed the same brightness. Although regions B and D 
had a similar brightness in August 2004, region D was ~ 0.9 magnitudes fainter (AK=2 vs. 
AK=0.7) than region B in July 2007. The knot that is pointed out in Figure 3 is a feature 
that reappears in the same place after a period of not being visible. It is faintly visible in 
our July 1999 observation, but it is much more prominent just three months later. It is not 
visible in 2003, but it is again prominently visible in May 2004. By July 2007, it is only 
faintly visible. 

We note that in the course of our observations, neither the position nor the morphology 
of these knots, as well as the features around them, ap pear to change. Rather, only their 



apparent brightness changes. Similarly, iLamplandl (119261 ) noticed that morphological features 



in Bubble's Variable Nebula would disappear only to once again become visible in the same 
place. This led to his conclusion that the most likely explanation for the observed changes 
in morphology is changes in the illumination of the nebula. Since we observed very similar 
phenomena in IRAS 18148-0440, we propose that the cause of the morphological variability 
of the nebula are changes in the illumination of the nebula, possibly by opaque clouds very 
close to the source. As noted above, morphological changes are noticeable within two months, 
with significant changes in the illumination of the cloud are observed to occur within a year 
(e.g. between July 2003 and May 2004). If w e assume that opaque clouds in Keplerian orbits 



around a 0.35 Mq ( iFuUer fc Woottenll2000l ) star are the cause of the change of illumination 
of the nebula, then the time scale of the observed variability corresponds the orbital period 
of the clouds. In this case, a period of a year suggests that the distance from the central star 
to these clouds is of the order of ~ 1 AU. However, images taken several years apart show 
a similar pattern in the differenced images (e.g. Aug 2004 to Jul 2007), suggesting that the 
time scale for significant changes in the illumination of the cloud may be several years, in 
which case the obscuring clouds would be farther from the central star. 



Since this is one of a very few known infrared variable nebula, it is difficult to estimate 
how common such nebulae are. Since there are only a few optically visible variable nebulae 



- 8 - 



that are illuminated by T Tauri stars, and since there are an order of magnitude fewer Class 
I YSOs than T Tauri stars, infrared variable nebula may be very rare. However, Class I 
YSOs are much more deeply embedded than T Tau ri stars, and commonly have a reflection 
nebula visible in the near-IR (jConnelley et al.ll2007l ). As such, future observations may show 
that there are variable near-IR nebulae around many other Class I YSOs. 



3.4. Ho Emission 



Figure 6 presents 2' wide images of IRAS 18148-0440 to illustrate the relationship 
between the near-IR reflection nebula, the H2 knots, as well as the location of the source 
and the outflow cavity which are visible in the mid-IR. The bright H2 shock to the west, as 
well as the H2 knots farthest east, are visible on account o f thermal dust ern ission in Spitzer 
IRAC images at 3.6 jjm, 4.5 /um, 5.8 /zm, and 8.0 fim (lEvans et al.l 120031 ) . The infrared 
emission corresponding to the H2 knots are found 65" to the west and 118" to the east 
of the so urce, the loc a tion o f whi ch is clearest i n the 8.0 fim Spitzer image. However, CO 
maps by iFuUer et al. Jl995h and iTafalla et al.l (120001 ) show that the CO outflow, which is 
aligned with the axis of the H2 knots, extends about as far east of the source as to the west. 
Since the CO outflow extends only as far as the bright H2 shock to the west, t he eastern H^ 



knots are farther east from the source than the CO emission contours shown in Tafalla et al. 



( I2OOOI ). We note that the western knots appear brighter, likely because the western side of 
the nebula corresponds with the blue shifted lobe of the CO outflow. The western knots 
are thus seen through less extinction than the H2 knots from the more deeply embedded red 
shifted lobe of the outflow. None of the H2 knots correspond with conspicuous knots in the 
bright reflection nebula, and no H2 knot is found within 20" of the source. 



3.5. Nebula Morphology 



Fuller et al.l (119951 ) showed that the near-IR reflection nebula is coincident and aligned 
with the blue shifted lobe of the CO outflow from the source. The near-IR image of the nebula 
shows the half of a bipolar cavity that is closer to us, where flux from the central source is 
scattered by the inner walls of the bipolar cavity. It is quite common for only one half of 



a (pr esumably) bipolar cavity to be visible in an infrared reflection nebula (jConnelley et al. 



20071 ). Most, if not all, of the flux in the 8 /im image is scattered light since the FWHM of 
the bright source in the center of the nebula is 4"22 whereas the FWHM of a nearby star 
is 1"91. The centroid of the 24 ^m source observed by Spitzer is 1"1 to the west of the 
8 /Ltm centroid. At 24 /im, this source is observed to be unresolved since the FWHM of the 
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24 /iin source is consistent with the FWHM of nearby field stars. This is hkely to be the 
true location of the protostar, and this location is shown in Figure 3. 

The dark lane visible ~2'.'5 (~560 AU) west of the bright 8 /im source appears to be 
coincident with the sharp eastern edge of the near-IR refiection nebula. This dark lane 
appears to be similar to a disk shadow being projected into the cloud (Hodapp et al. 2004, 
Pontoppidan et al. 2005). If this dark lane was a disk shadow, then the brightest part of 
the nebula should be on the closer (blue shifted) side of the bipolar cavity since forward 
scattering of light by dust is more efficient than back scattering, and light scattered on the 
closer side of the nebula suffers less extinction. However, the bright 8 /xm source is on the 
farther (red shifted) side of the nebula, which is difficult to explain within the conventional 
model of a single source with a disk illuminating a bipolar cavity. It is possible that the 
illuminating source of the nebula is totally obscured and the dark lane is a shadow of the 
disk being projected into the cloud, which would account for the dark lane being prominent 
from 1.6 fim to 8 fim. In this case, the 8 /xm source may be scattered light from a nearby 
and less obscured companion. 



4. Summary 

We have presented K-band images of IRAS 18148-0440 in L483 spanning 17 years. This 
ClassO/I YSO is one of the first large nebulae that are only visible in the near infrared that 
has been found to be variable. Morphological changes are noticeable within two months, and 
features are observed to disappear then reappear in the same location. These changes are 
consistent with a changing illumination of the nebula, which may be caused by opaque clouds 
within ~ 1 AU the protostar. As such, this nebula is an infrared analogue to well known 
optically visible variable nebulae such as Bubble's Variable Nebula. The overall brightness 
of the nebula appears to be slowly declining, with observed variations of a magnitude within 
two months and two magnitudes within a year. Since these changes happen on a time scale 
a few times greater than the light crossing time of the nebula, the intrinsic changes in the 
illumination of the nebula must be rapid. The H2 shocks are visible in archival Spitzer 
images, and these shocks can be found nearly twice as far to the east of the source as to 
the west. The SED of this sou rce is consistent with it be ing a Class I YSO on account of 



its near-IR fiux. While Cep A (IHodapp fc Bressert II2OO8I ) and L483 host the first infrared 



variable nebulae identified to date, these are likely to be the first of many to be found. 
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Fig. 1. — The SE P of L483, combin i ng da ta from 2MASS, Spitzer, IRAS and an 850 /im 
observation from [Fuller &: Wootteru (120001 ) . There is significant emission at wavelengths 
shortward of 10 /im, and the SED is not well fit by a single temperature blackbody, which 
is uncharacteristic of a Class YSO. The dotted line is a blackbody with a temperature of 
45 K, which was temperature of the best-fit blackbody function to the IRAS data. We did 
not fit a separate blackbody function to the flux values from 1.2 //m to 8.0 /xm since fluxes 
in that wavelength range are strongly affected by extinction and scattered light. 
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Fig. 2. — The variability of the integrated flux of the L483 infrared nebula at K-band. The 
plus, diamond, triangle, and 'x' symbols represent photometry of the nebula relative to four 
different field stars. The square is the K-band magnitude measured by 2MASS. The dotted 
line passes through the average value at each epoch. 



-14- 



B C knot 




D 



Fig. 3. — The master image with labels for four regions as well as a knot. Changes in the 
brightness of these regions are discussed in section 3.2. The location of the 24 yum source 
observed by Spitzer is marked with the star symbol. 
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Fig. 4. — 30" by 20" views of L483 at K or K' at 8 epochs taken over a span of 17 years. 
The grayscale was chosen to best show the morphology of the nebula, and does not represent 
the actual surface brightness of the nebula. The left column shows the images of the nebula 
whereas the right column shows the difference between the image on the left and the master 
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Fig. 4. — Figure 3 continued. 
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Fig. 5. — K-band photometry of four regions within the nebula. Data from October 1990 
is to the left of the plot for clarity. Generally, different regions of the nebula brighten and 
fade together, occasionally by different amounts. However, there are notable exceptions, 
particularly between 2003 and 2004. The image taken in July 2002 is not deep enough to 
allow us to measure the photometry in different regions of the nebula. 




Fig. 6. — 120" by 60" views of L483 with North up and East to the left. The top frame 
is a K-band image from July 1999. The source lies behind the small, faint nebula between 
the main nebula and the fainter nebula to the East. Note the dearth of field stars within 
~ 45" of the source. The middle frame is our H2 image with the K-band image from the 
same night subtracted out. The bright shock to the west is also visible in the K-band image. 
The H2 shocks extend farther to the east than the west from the location of the source. The 
bottom frame is a Spitzer IRAC 8.0 //m image. Both lobes of its outflow cavity are now 
clearly visible. The bright source near the center of the nebula has over twice the FWHM of 
field stars, and thus is likley scattered light from the central star, which remains obscured. 
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Table 1. Observations 



Date (UT) 


Observer 


Filter 


FWHM (") 


Instrument 


Oct. 6, 1990 


Hodapp 


K' 


0.92 


NICMOS3 


June 27, 1993 


Puller 


K 


2.68 


IRIM 


Apr. 15, 1999 


2MASS 


K 


2.31 


2MASS 


July 25, 1999 


Hodapp 


K 


0.52 


QUIRC 


July 25, 1999 


Hodapp 


Ha v=l-0 S{1) 


0.64 


QUIRC 


July 26, 1999 


Hodapp 


H 


0.60 


QUIRC 


Sep. 21, 1999 


Hodapp 


K' 


0.55 


QUIRC 


July 3, 2002 


Hodapp 


K' 


1.05 


QUIRC 


May 14, 2003 


Connelley 


K 


1.15 


QUIRC 


July 18, 2003 


ConncUcy 


K 


1.38 


QUIRC 


May 28, 2004 


Connelley 


K 


0.70 


SpeX 


Aug. 1, 2004 


Connelley 


K 


1.26 


QUIRC 


July 13, 2007 


Connelley 


K 


1.02 


SpeX 



